Electron fractionalization in two-dimensional graphenelike structures 
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Electron fractionalization is intimately related to topology. In one-dimensional systems, fraction- 
ally charged states exist at domain walls between degenerate vacua. In two-dimensional systems, 
fractionalization exists in quantum Hall fluids, where time-reversal symmetry is broken by a large 
external magnetic field. Recently, there has been a tremendous effort in the search for examples of 
fractionalization in two-dimensional systems with time-reversal symmetry. In this letter, we show 
that fractionally charged topological excitations exist on graphenelike structures, where quasiparti- 
cles are described by two flavors of Dirac fermions and time-reversal symmetry is respected. The 
topological zero-modes are mathematically similar to fractional vortices in p-wave superconductors. 
They correspond to a twist in the phase in the mass of the Dirac fermions, akin to cosmic strings 
in particle physics. 
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In low dimensional systems, the excitation spectrum 
sometimes contains quasiparticles with fractionalized 
quantum numbers. A famous example of fractionaliza- 
tion was obtained in one-dimension (ID) by Jackiw and 
Rebbi [l[ and by Su, Schrieffer and Hegger Q. They 
showed the existence of charge e/2 states, with poly- 
acetelene as a physical realization of such phenomena. 
In these systems, a charge density wave develops and the 
ground state is two-fold degenerate. The fractionalized 
states correspond to mid-gap or zero-mode solutions that 
are sustained at the domain wall (a soliton) interpolating 
between the two-degenerate vacua. 

The fractional quantum Hall effect provides an exam- 
ple of fractionalization in two-dimensions (2D). Not only 
do the Laughlin quasiparticles have fractional charge 
but they also have fractional (anyon) statistics 
Time-reversal symmetry (TRS) is broken due to the 
strong magnetic field, leaving as an outstanding problem 
the search for systems where fractionalization is realized 
without the breaking of TRS. The motivation for such a 
quest stems from speculations that fractionalization may 
play a role in the mechanism for high-temperature su- 
perconductivity @, 0, @|- Progress has been made on 
finding model systems, such as dimer models, in which 
monomers defects act as fractionalized (and deconfined, 
in the case of the triangular lattice) excitations. 

In this letter, we present a mechanism to fractionalize 
the electron in graphenelike systems that leaves TRS un- 
broken. The excitation spectrum of honeycomb lattices, 
which have been known theoretically for a few decades to 
be described by Dirac fermions is now the subject 

of many recent studies since single and few atomic-layer 
graphite samples have been realized experimentally [11 1 . 
Quantum number fractionalization is intimately related 
to topology and here we find that a twist or a vortex in 
an order parameter for a mass gap gives rise to a single 
mid-gap state at zero energy. Such twist in the mass of 



the Dirac fermions in graphenelike structure is the anal- 
ogous in 2+1 space-time dimensions of a cosmic string in 
3+1-dimensions [l2T |. 

The zero-modes we find are, in their mathematical 
structure, similar to those found in p-wave superconduc- 
tors by Read and Green [13[ and to those found by Jackiw 
and Rossi [14| when demonstrating topological excita- 
tions and suggesting an index theorem in 2D (see also 



Ref.U 



In Refs. 



14 , and |15|, a twist in the phase 



of a superconducting order parameter and a charge 2e 
Higgs boson, respectively, were considered and thus the 
electric charge assigned to the zero-mode is not a good 
fractional quantum number. Instead, in the systems we 
study, the electronic charge is conserved and therefore we 
can show that it is fractionalized by properly tallying it. 

We find only one normalizable state for each vortex, de- 
spite the two flavors of Dirac fermions in honeycomb lat- 
tices. The fact that there is one and not two zero modes 
is essential for fractionalization: doubling the number of 
zero modes due to two flavors of Dirac fermions, in ad- 
dition to the spin degeneracy, would lead to excitations 
with the same quantum numbers as ordinary electrons. 
We find excitations with charge Q = e/2 in the spinless 
case, and charge Q = e and spin S = or charge Q = 
and spin S — 1/2 in the case with spin. Hence, these 
fractionalized quantum numbers are similar to those in 
polyacetelene, but the counting is different in ID and 2D, 
and as we will discuss, there is a C/(l) symmetry underly- 
ing the fractionalization in the 2D case which would not 
lead to fractionalization in ID (in which a Z 2 symmetry 
is underlying the effect). 

The mechanism for fractionalization can be portrayed 
in a simple form by considering spinless electrons hop- 
ping on a honeycomb lattice with textured tight-binding 
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FIG. 1: (color online). The honeycomb lattice A with lattice 
spacing a and its two triangular sub-lattice A A (Red) and A fl 
(Blue). The generators of A A are ax and a 2 . The three vec- 
tors Sj connect any site from A A to ist three nearest-neighbor 
sites belonging to A s . The Kekule distortion is a modulation 
of the nearest-neighbor hopping amplitude that is indicated 
by representing nearest-neighbor bonds of the honeycomb lat- 
tice in black (grey) if the hopping amplitude is large (small). 



hopping amplitudes and described by the Hamiltonian 



H 



reA, i=l 



H.c. 



(1) 



Here, s i with i — 1,2,3 connects any site r belong- 
ing to the triangular sublattice to its three nearest- 
neighbors belonging to the triangular sublattice A B of the 
honeycomb lattice as is depicted in Fig. [T] The fermionic 
annihilation operators a r and b r act on A^ and A B , re- 
spectively, and so do their adjoints. Graphene is often 
described by Hamiltonian |1} in the single-particle ap- 
proximation, neglecting the spin of the electron, and as- 
suming 5t ri = 0. We are going to show that the small 
variations of the hopping strength, St r i , over the uniform 
hopping t provide the background on which fractionally 
charged states can be constructed. We then discuss how 
such 8t r i can arise from a local order parameter that 
decouples electron-electron interactions. 

When St r i = 0, Hamiltonian ([1]) can be diagonal- 

ized in momentum space, H = J2k^k a k^k + H.c, 



* h = -*£;=!■ 



The single-particle spectrum e k — 



± thus contains two (zero-energy) Dirac points at 



the zone boundaries K 



± 



± 



(- 



, ) . After lineariza- 



tion, k — K ± +p, one obtains a spectrum containing two 



chiral flavors ±. 



"p a p,± v P 



H.c, 



y P,± 



ivp(p x Hpy), with the Dirac conelike structure £±(p) 



±v F \p\ for the energy dispersion. (Hereafter, 



1.) 



We shall focus on backgrounds that yield a chiral mix- 
ing between the ± species. A Kekule texture, depicted 
in Fig.[TJ provides such a mixing [l6j |. The modulation 



St ri = A(r) e iK +' Si e lG r /3 + c.c, 



(2) 



with wave vector G := K + — K_ couples the Dirac 
points at K± as is depicted in Fig. O Here, we al- 
low for spatial fluctuations (on length scales much longer 




FIG. 2: (color online). The first Brillouin zone of the tri- 
angular lattice is shown together with the Dirac points K± 
at the zone boundary, the reciprocal vector G that connects 
them, and the massive relativistic dispersion centered about 
K± that opens due to a Kekule distortion with A(r) = A . 



than the lattice spacing a) of the complex-valued order 
parameter A(r). The phase of A(r) controls the or- 
dered hopping texture. To leading order in a gradient 
expansion, Hamiltonian ([1]) subjected to the Kekule tex- 
ture © is given by H = J d 2 r ^(r) 1C D (r) *(r) with 
*t( r ) = ( u t( r ) u t( r ) „t( r ) vl( r ) ) and 



( ~2id z A(r) ^ 

~2id g A(r) 

A(r) 2id z 

\ A(r) 2id B / 



(3) 



We are using the notation z = x + iy, d z = (d x — id y )/2, 
with an overline to denote complex conjugation. With- 
out the Kekule texture ([2]) leading to A(r), the Dirac 
kernel (|3|) is the single-particle relativistic massless Dirac 
Hamiltonian in 2+1 space-time. With the Kekule tex- 
ture ([2]) A(r) = A , the dispersion takes the simple form 
£±{p) = ±\/|p| 2 + l^ol 2 ' i- e -' a single-particle mass gap 
|A | has opened. The Dirac kernel © is TRS. TRS orig- 
inates in the tight-binding hopping elements being real. 
Moreover, the transformation law under a r — > —a r and 
b r+s — > +b r+s of the single-particle tight-binding ([1]) 
ensures that any positive energy eigenstate of the Dirac 
kernel ([3]) can be matched to a negative energy eigen- 
state, while only zero-modes can be left unmatched. We 
shall call this property sublattice symmetry (SLS) (l7j . 

The order parameter A(r) = A can be complex val- 
ued, but he spectral mass gap is real. This suggests that 
the phase of A is redundant. In fact, it can be removed 
from Eq. ([3]) with a chiral transformation that rotates the 
phases of the ± species by opposite angles. This is not 
true anymore if the phase of the order parameter A(r) 
varies in space, and, in particular, if it contains vortices. 
The latter situation leads to fractionalization. 

We assume that 



A(r) = A (r)e 



i(a+n6) 



(4) 



where A (r) > 0, n £ Z for a single- valued order parame- 
ter, and we are using the polar coordinates z = rexp(i^). 
We are seeking eigenstates of the Dirac kernel © with 
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vanishing energy that are normalizable, i.e., 

(a r -zr- 1 a e ) Ua (r)+ le l9 A(r) v a (r) = 0, (5a) 
ie- l6 A(r) Uo (r) - (d r + ir-%) v a (r) = 0, (5b) 

holds on sublattice A A while two more equations ob- 
tained from Eq. |5j) with the substitutions u a — > u b , 



u fc , and 



-0 must also hold on sublattice A, 



The same equations were studied in Refs. [H, [3 from 
a different perspective since the textured gap in these 
works is associated to a superconducting order parame- 
ter. Hence, electric charge is not a conserved quantum 
number in Refs. 13l. 1 141] . There are \n\ independent nor- 



malizable zero modes, which are either supported on sub- 
lattice A a when n < -loron sublattice Ab when n > 1. 

We assume that n = —1 in which case the single- valued 
and normalizable wave functions for the zero mode is 



u a (r,6) 



*(f + f) 



f;dr' A (r') 



20F v^T^ 



-2/J-dr'A (r') 



(6) 



Its support is on sublattice A A . The zero mode when 
n = 1 is obtained from Eq. (|6]) with the substitutions 

u a( r >®) ~ * v b( r >®) an( i v a( r '®) ~ * u b( r '®)- When n = 1, 
the support of the zero mode is on sublattice A B . The 
wave function ([6]) decays exponentially fast away from 
the core of the vortex ((4| at the origin of the complex 
plane. Its localization length is set by the gap value A 
reached when Aq(t') saturates sufficiently far from the 
origin. If the Kekule texture ^ supports a pair of vor- 
tices a large distance R apart, then the Dirac kernel ([3]) 
has two eigenstates whose energy eigenvalues are expo- 
nentially small. 

To obtain the charge bound to a vortex, one has to 
study 5v(r,e) = V\ n \ =1 (r,e) - ^ n | =0 (r,e) where v n (r,e) 
is the local density of states (LDOS) of the Dirac ker- 
nel ([3]) in the presence of the mass twist Q. Because 
of the SLS, to any negative eigenstate of the Dirac ker- 
nel, tp_ E (r), there corresponds a positive energy state, 
ijj +£ (r), related to ijj_ £ (r) by a unitary transformation. 

Hence, the LDOS v{r,e) = £ E , ^,(r)^ £ ,(r)5(e - e') is 
symmetric with respect to zero energy and negative and 
positive energy eigenstates contribute equally to the local 
density v(r) = lim £ _,<x) J^ e ds'v(r, e'). In the presence 
of the zero mode tpo( r ) this spectral symmetry together 
with the conservation of the total number of states imply 



d 2 r 




5v(r,e)de + \Mr)\ 2 = °- ( 7 ) 



But the single zero mode tpo( r ) ^ s normalized to one and 



/ d 2 r [ 8v{r,e)de = -1/2. 

J J — OO 



(8) 



According to Eq. ([8]), the valence band has a deficit of 
half of a state, as does the conduction band. Thus the 
difference in net charge between a fully occupied valence 
band with and without the vortex is — e/2. The total 
charge of a closed system must still be an integer. This 
can be understood as follows. Vortices must be present 
in pairs (for example, if a finite system is to satisfy peri- 
odic boundary conditions or for energetic reasons in the 
thermodynamic limit), and while the charge around each 
vortex is half that of an electron, the fact that the vor- 
tices appear in pairs ensures that the total charge is an 
integer multiple of e. The mechanism for fractionaliza- 
tion is similar to that in ID [l|| , but there the fractional 
charge stands in domain walls between two degenerate 
vacua, which must appear in pairs (again if the system is 
to satisfy periodic boundary conditions). In ID, however, 
one needs a discrete symmetry to fractionalized charge. 
For example, for charge e/2 states one needs a Z 2 sym- 
metry, which can be understood generically in terms of 
an accumulation of charge along phase twists that have a 
natural interpretation within a bosonization scheme [l^ ] . 
These arguments do not extend to 2D, which is our con- 
cern in this letter, as there is no simple connection be- 
tween phase twist and charge as in the ID bosonization 
scheme. The fact that the fractional states must occur in 
pair, nonetheless, follows from the common requirement 
for the Z7(l) mass vortices in 2D and the Z 2 domain walls 
that these defects occur in pairs to satisfy appropriate 
boundary conditions. 

In the case of electrons with spin, two independent 
zero modes that carry the spin quantum number are ex- 
pected. Hence, there is, in total, a full electronic state 
missing in the valence band for each vortex. Conse- 
quently, unoccupied, singly, and doubly occupied zero 
energy states correspond to the charge and spin quan- 
tum numbers (Q = -e, 5 = 0), (Q = 0, S = 1/2), and 
(Q = +e, 5 = 0), respectively. Notice that, had we ob- 
tained two normalizable solutions for a vortex or an anti- 
vortex as opposed of one as in Eq. (|6|), all the quantum 
numbers would be "doubled" and simply coincide with 
those of ordinary electrons. 

These arguments for fractionalization require existence 
of the mid-gap state and the SLS of the spectrum. SLS 
is broken when a next-nearest neighbor hopping t' is in- 
troduced. A t' brings an overall shift and a quadratic 
term into the diagonal sector of the Hamiltonian, 87i' = 
(3t' + ^-j-p 2 ) 1, up to second order in p. This pertur- 
bation alters the energy spectrum, but Dirac points (in 
the absence of the background) remain when the systems 
is at half-filling. The textured background still opens a 
gap at the Dirac points and there is a mid-gap state if 
there is a mass twist. However, the mid-gap state is not 
exactly at the center between the valence and conduction 
bands, being shifted by Se^ = ^^-Ag to first order t'. 
As long as this shift is small compared to the gap, the 
mid-gap solution is robust against breaking of the SLS. 
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The fractionalized quantum number persists along with 
the single bound state, but the argument for charge e/2 
is not valid. Starting from zero t' , turning it adiabati- 
cally, one can argue that fractionalization should remain 
but take irrational values as in id dim!!!. 

We now turn to a mechanism for generating sponta- 
neously the Kekule distortions from repulsive electron- 
electron interactions. We assume the nearest- neighbor 
interaction H' := X^eA A a r & t-+ Sj a r 

with V the interaction strength. Evidently, this inter- 
action preserves TRS, the point group symmetry of the 
honeycomb lattice, and SLS (up to the total number 
operator). After linearization around the unperturbed 
Dirac points, a mean-field decoupling of this interaction 
with respect to the order parameter 

-^A=(bl(r + s j )a + (r)) = (al(r) b + (r + Sj )) (9) 

coincides with the case of the Kekule distortions. If the 
system is at half-filling and zero temperature, all neg- 
ative energy states are occupied and the self-consistent 
^ cxp(i P . aj .) follows ^ The gum can 



equation 1 



2|AJ 



be replaced by an integral up to an appropriate momen- 
tum cutoff, Aa ~ 33/i 1 ncre chosen so as to match 
the total number of states in the effective theory to 
that in the microscopic theory. As this integral con- 
verges even without the |A| 2 term, neglecting it yields 
1 < 3v f|^ A ° [Jo( A «) H -i(^) - J_i(Aa) H (Aa)]. Here, 
J u (z) is a Bessel function and H v (z) is a Struve function. 
Insertion of the numerical value for Aa then yields the 



lower bound V„ > 



1*1 



1.906|i|. Self-consistency de- 



0.52465 _ 

mands a repulsive interaction. (The very same mean-field 
theory also renormalizes the uniform hopping: t — ► t+St, 
St ~ 0.1V.) We note that the condition for the Kekule 
pattern to form could potentially be engineered in a gas 
of dipolar fcrmionic atoms trapped in an optical honey- 
comb lattice, where the ratio V/t could be tuned. 

At last, we need a vortex on top of the Kekule distor- 
tion {2} for fractionalization to happen. Here, we face 
a difficulty in that the chiral symmetry of Eq. (|3|) is an 
artifact of the linearization. Integration of the fermions 
on the lattice yields an effective action for the phase of 
A with a discrete symmetry Z 3 . Vortices are then con- 
fined at sufficiently long distances due to a pinning po- 
tential for the phase of A. Fractionalization can only 
be observed on length scales much larger than A yet 
much smaller than the confining length scale, a func- 
tion of A /i. This difficulty can be overcome by a small 
breaking of the point-group symmetry of the honeycomb 
lattice by assuming anisotropic hopping between nearest- 
neighbor sites. By continuity, the only effect of such a 
small anisotropy is to move the Dirac points away from 
the boundaries of the first Brillouin Zone, as is confirmed 
by an exact computation of the spectrum in Ref . [23|. The 



Kekule wave vector G in Eq. {2]| is then not anymore 
commensurate to the reciprocal lattice, integration over 
the fermions does not produce a pinning potential for the 
phase of Kekule modulation, and vortices can prolifer- 
ate as a result of temperature-induced fluctuations (their 
bare logarithmic interaction is screened for temperatures 
above the Kosterlitz-Thouless transition temperature). 

In summary we presented a mechanism in graphene- 
like 2D condensed matter systems that realizes charge 
fractionalization with TRS. The mechanism is based on 
an effective low-energy Hamiltonian of the Dirac type 
with a textured mass with phase twists due to vortices. 
When the vortices are far apart, each vortex carry a zero- 
mode responsible for local charge fractionalization. This 
mechanism is a condensed matter 2+1 space-time ana- 
logue to 3+1 space-time cosmic strings. 
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ported by the NSF grant DMR-0305482 (C.-Y. H and 
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